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We give a detailed account of the free field spectrum and the Newtonian limit of the 
linearized "massive" (Pauli-Fierz), "topologically massive" (Einstein-Hilbert-Chern- 
Simons) gravity in 2 + 1 dimensions about a Minkowski spacetime. For a certain ratio 
of the parameters, the linearized free theory is Jordan-diagonalizable and reduces to 
\ a degenerate "Pais-Uhlenbeck" oscillator which, despite being a higher derivative 

' theory, is ghost-free. 
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^ . As is well known, Einstein's General Relativity (GR) in 3+1 dimensions is nonrenormal- 

Q \ izable and breaks down at "high" energies. If GR is considered as a low energy effective 
OO ■ limit of some unknown fundamental theory, then it will receive higher curvature corrections, 

^ ■ such as al^ + PRf^^, etc. which give a much well-behaved gravity theory in the UV region, 
O ■ but unfortunately ghosts are not decoupled and unitarity is lost. In fact, in the Einstein 

' plus quadratic curvature model whose Newtonian potential reads 
^■ 

V{r) = Gm{ - i - t'- + 

I r 6 r 6 r 

one can see the nice UV behavior (V(0) = 0) and the presence of the ghost related to the 
last repulsive term (see [1] and the early references therein). It is, therefore, of extreme im- 
portance to get a higher derivative gravity theory without ghost problems. 2+1 dimensional 
gravity, which is quite easy compared to GR, is a nice laboratory to study various quantum 
gravity issues. Admittedly, the lesson one learns for 3+1 dimensions should always be taken 
with a grain of salt. Nevertheless, valuable insight on higher derivative gravity theories, as 
well as massive gravity theories, can be gained with the study of 2+1 dimensional gravity 
models. In fact, an analogous approach in quantum field theory has borne much fruit. In 
this paper, we shall find a particular ghost-free higher derivative gravity model, which has 
various "mass" terms. 

In 2 + 1 dimensions, in addition to the usual gauge non- invariant Proca/Pauli-Fierz mass 
terms, one can add a gauge invariant "mass" which is the Chern-Simons term [2] to both 
(the abelian and the non-abelian) spin-1 and spin-2 models. [To be more precise, the Chern- 
Simons term acts like a mass term only at the quadratic level.] There are major differences 
between the vector and the tensor cases, one of which is related to the gauge invariance: In 
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the former, even the Proca mass can be made gauge invariant by introducing a Higgs field 
[see e.g. [3] and the references therein]; but in the latter, it is simply not known how to get 
the Pauli-Fierz mass from a gauge (diffeomorphism) invariant action. 

The van Dam-Veltman-Zakharov ([4, 5]) discontinuity between the massive and massless 
gravity models registers as another difference between the vector and tensor theories. [Recall 
that the vDVZ discontinuity states that around the Minkowski background, the massless 
limit of the massive gravity (Einstein- Pauli-Fierz theory) is not the massless gravity (Einstein 
theory).] The existence of such a discontinuity [15] and perhaps the isolation of massless 
gravity from the massive (albeit arbitrarily small) ones, is of extreme importance as far as 
the concept of mass in gravity is concerned. The discontinuity also has much to do with the 
difference in the degrees of freedom between the massive and massless models. 

In a previous Letter [7], the 3-term model (Einstein-Hilbert, Chern-Simons and Pauli- 
Fierz) was studied in the context of various discontinuities and parameter ranges. The 
degrees of freedom were identified but the diagonalization of the model (or the explicit 
form of the free field Lagrangian) was not carried out. As we shall see in this current 
work, a detailed analysis of the free theory is not redundant at all, since the theory exhibits 
rather interesting features: For example, as opposed to being a collection of simple harmonic 
oscillators, we shall get higher order degenerate (i.e. equal frequency, mass) Pais-Uhlenbeck 
oscillators [8] and the model (for a certain ratio of the Chern-Simons and the Pauli-Fierz 
parameters) can be brought into a Jordan-diagonal form at most. This is rather unexpected 
since, without the Pauli-Fierz mass term, the free field limit of the model is not a higher 
derivative model. Thus we have a very interesting situation here: the higher derivative model 
is ghost-free! In fact, in the recent literature, one can find examples of similar ghost-free 
higher derivative models: Mannheim-Davidson [9], Smilga [10, 11] and Hawking- Hertog [12] 
all deal with the one- dimensional degenerate Pais-Uhlenbeck oscillator as a model of higher 
derivative ghost-free theories. 

There are a number of reasons why one would be interested in higher derivative models: 
First of all, they are much better behaved in the UV region (see e.g. Stelle's [13] work 
in higher curvature gravity) and secondly, even if one does not have higher order terms at 
the classical level, quantum corrections usually introduce such terms. Ironically, in most of 
the renormalizable theories, during the process of regularization and renormalization, one 
introduces ghosts only to remove them at the end! Thus a higher derivative model without 
ghosts, as the one considered here, would be of some interest. 

The layout of the paper is as follows: In section 2, we introduce the 3-term model and work 
out the free field Lagrangian. In section 3, we suppress the spatial dependencies to basically 
look at the model as a collection of oscillators. In the first appendix, we give a basic review 
of the Pais-Uhlenbeck oscillator and the second appendix contains the Newtonian limit of 
the 3-term model. 

II. THE MODEL 

Since the signs and the values (or the ratios) of the involved parameters are of extreme 
importance, here we carefully define and work out the linear version of the 3-term gravi- 
tational action which is the sum of Einstein-Hilbert, (third derivative order) Chern-Simons 



3 



and standard Pauli-Fierz mass, terms 

^ (1) 

only at quadratic order in deviations h^u = g^^ — Vtiu about the flat Minkowski background 
ry^i/ with the usual definition h = rj^^^h^^. Our signature is (—,+,+), we use eoi2 = 1 
throughout and, in what follows, all operations are carried out with respect to the flat 
background ?7^j,. The model under investigation is the most general such model [7] since 
the sign of fi is arbitrary, that of is free and the presence of a allows for choosing the 
Einstein term's sign (a = +1 in the usual case, a = —1 for topologically massive gravity 
(TMG) [2] and a = for CS-Pauh-Fierz theory). Note that — > oo represents massive 
gravity with 2 excitations and m = is TMG with a single mode [2] , whereas pure Einstein 
theory, 1/fi = = 0, has no excitations in D = 3. Without working the details, a cursory 
look at the model does not really reveal what kinds of "oscillators" we have for the free field 
limit. In fact, as will be seen shortly, there are various discontinuities as far as the limits of 
the parameters are concerned. 

We shall be interested only with the free field limit for now. [Of course, one eventually 
has to deal with the interacting theory, but, unfortunately, there are infinitely many terms 
in the interaction Lagrangian rendering the perturbation theory intractable.] The individual 
terms that make up the action (1) can each be expanded at quadratic order in h^i, as follows: 
The Einstein-Hilbert piece reads 

lE = aJ £x V^R =-^1 d'x V + 0{h') , (2) 

where 

= I + d^d'h'^ - ah'"' - d'^d'h) - \ri"' {do^dph'"'^ - nh) . (3) 

The Chern- Simons bit yields 

^cs^-^J d^xe^'-'T^, (^0^r%, + h\f3T^,,^ "~^/ d'xe^apGrd'^hf'^ + Oih'). 

(4) 

Finally the Pauli-Fierz mass term, up to 0{h^), gives 

It was shown in [2] that the sign of the Einstein term in pure TMG (m^ = 0) must 
be a = — 1, opposite to that in the usual Einstein gravity, in order for the energy to be 
positive, independent of the sign of On the other hand, the usual massive spin-2 system 
does have excitations and this forces both the relative and overall signs of the Einstein and 
Pauli-Fierz mass terms to be in the usual way to avoid the presence of ghosts and tachyons. 
This indicates an inevitable dilemma in the choice of the sign of a in the TMG and massive 
spin-2 cases. 

To analyze the full system given in (1) for generic values of (//, m) and its dependence on 
the sign of a, we first decompose /i^i^ as 

hij = {Sij+didj)(f)-didjX+{eikdkdj+ejkdkdi)^ , hoi = -tijdjij+diNL , hoo = N, (6) 
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where di = di/ a/— 9|. The decomposition (6) yields the following for the Einstein (2), 
Chern-Simons (4) and Pauli-Fierz (5) components of (1), respectively: 



J d'x |(-V'r7 + O[V'(7V-27Vi)+x + n0]} , 



IpF = J d'x {NLV'NL + vy^v + e-<l>X + N((t> + x)} 



(7) 
(8) 
(9) 



When all these are put together, the terms proportional to can be collected altogether to 
obtain 

I ^^Jd'x |a((/)x-20V^7V^+(-V2r7 + O') + ^(-V^?7 + O[-2V^^i + X + □</>] 



aVV - -V\V^r] - - + x) 
A* 



(10) 



This looks hke a highly complicated system. To simplify it, the next thing to do is to 
use integrations by parts, whenever necessary, to eliminate those variables from (10) which 

can be treated as 'Lagrange multipliers' and to use the corresponding 'constraint equations' 
thus obtained so that only dynamical degrees of freedom are left. When N , r] and Nl are 
eliminated in this fashion, one finds 



I= - I d'x 



a(j)n(f) - m^(l)^ + 2XD(f) + a/i^X^ - m^UX - 0™ (M - - ^^C^ 



, (11) 



where wc have used A = — a0 + m^:^(0 + x) for simplicity. The remaining three dynamical 
components of h^i, yield 



50: DA + an0 - = , 

1 . . 1 . 



□0 - m^iZ—X + m^/i— ^ + ajjL^X = , 



5X : 

5^: /xA + ne = 0, 



which together give 



(□3 - aV^°^ + 2amV^n - mV)0 = , 
for the field 0, through which the remaining fields are determined as 



(12) 
(13) 
(14) 

(15) 



A 



—a + m 



and 



1 
□ 



1 



an — m n 



□ 



Equation (15) was obtained in [7] and [14], in which it was shown that for a = 1, the roots 
of the eigenvalue equation are complex unless yU^/m^ > 27/4. 

Introducing auxiliary variables t/^ = and Vt = 0^4> = Oip, the equation for (p (15) can 
be put in the form Dx — Ax, where 













1 





X = 




, A = 








1 




n 




4 2 

m ji 


— 2am^//^ 
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In passing, we should note that when ix^ /w? > 27/4, the eigenvalue equation has three 
distinct real roots, the explicit forms of which are rather cumbersome so that we refrain 
from presenting them here. However, the free field limit of the model can simply be thought 
of as being described by three uncoupled oscillators then. 

The degenerate case: 

Let us now consider the case when a = 1 and /i^/m^ = 27/4. Then det (pi — A) = 
(p — 3m^)^ (p — |m^) and by using the eigenvectors of A, one can form the modal matrix of 
A as 

16/9 1 
4mV3 Sm^ 1 
9m^ 6m? 



m 



which 'Jordan-diagonalizes' A as 



A = P-^AP = 



3m74 











1 

3m2 



By introducing z — P = [2^1 Z2 -^3]^, this system thus takes the form Dz — Az, which 
can be thought of as coming from an equivalent dynamical system whose action reads 



where specifically 



1 - 



3m^ 



□ + 



(fx 



zi{n- -m^ ] Zl + Z2{D - Zmyz2 



9m4 



I (p and Z2 



7 

-- + 



32 



9 27m2 



-□ - 



16 



81m' 



The first term represents the usual (ghost-free) massive real scalar field whose quanti- 
zation is well-known. The second term describes a degenerate (or purely quadratic) Pais- 
Uhlenbeck oscillator. In principle, just like any other higher derivative model, one expects 
this model to be contaminated with ghosts. However, as was shown in [9, 10, 12], this is not 
true: Purely quadratic models differ from those with quadratic plus lower derivative ones. 
This should not be surprising because, even the classical solutions of the purely quadratic 
theories have milder instabilities and physically acceptable parameter ranges. 

This, of course, is all at the level of the free-field theory. When interactions are introduced, 
the picture will necessarily change drastically. We are not aware of much work on this, 
save the work of Smilga [11] who has shown that once a quartic interaction is introduced, 
depending on the numerical values of the ratios of the involved parameters (such as the 
coupling constant to mass ratio), the model still has ghost-free, stable regions. This issue 
obviously deserves more attention. 



III. THE REDUCED MODELS 



In this section, we shall suppress the spatial part of the theory to see more transparently 
the degrees of freedom (the free oscillators) of the model. The advantage of this is that (10) 
simplifies a great deal without losing any of its degrees of freedom. Much of what we have 
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done above (in a fully relativistic way) can be understood more easily from the following 
reduced Lagrangian 

^ = ljd'x |a(0x + e') + ^e(x-0)-m'(f -0X + A^(0 + X))} • (16) 
Getting rid of the Lagrange multiplier N, we find the reduced form of the 3-term model: 

I^IJ |a(x^ + e) + Ux-m\e + X') } , (17) 

which leads to the following equations of motion: 

X + m\^ - e , C + m^C = . (18) 

Let us now check various limits: 
The pure TMG case: 

Setting the Pauli-Fierz mass to zero (m^ = 0) and eliminating ^, which behaves like a 
Lagrange multiplier, we have 

I=\jd'x{-a{x'-a'l^\')]. (19) 

It is clear that when a — —1, we have a ghost-free, massive, single excitation. This is the 
good old result of [2]. 

Einstein-Pauli-Fierz case: 

Setting /X — >• oo, we get 

I^ljd'x [a{x' + i') - m\e + X')]. (20) 

which describes two real scalar, massive degrees of freedom if a = +1. [This is what one 
expects from the onset, since massive spin-2 field has 2 degrees of freedom in 2+1 dimensions, 
that is as much as the massless spin-2 field in 3-1-1 dimensions.] 
No Einstein term: 

Setting a = yields a higher derivative model 



) , (21) 



which has to be diagonalized. However, as the eigenvalue equation (M^ + |yu|m^ = 0), where 
M is the mass of the excitations, constructed from the equations of motion, shows, this 
model has tachyonic solutions. This could be of some interest in Anti-de Sitter spacetime 
backgrounds which allows for certain negative mass squared solutions, but is not viable in 
flat backgrounds. 
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IV. CONCLUSIONS 

In this work, we have studied massive gravity (with two different types of mass: the 
Chern-Simons and the Pauli-Fierz) in 2+1 dimensional gravity. The theory has interesting 
mass discontinuities: It becomes a higher order (degenerate) Pais-Uhlenbeck oscillator if the 
ratio of the Chern-Simons and the Pauli-Fierz masses is tuned to 3v/3/2. The full theory 
(i.e. the Einstein-Chern-Simons-Pauli-Fierz model) generically has, at the linearized level, 
3 distinct massive degrees of freedom. The free field Lagrangian is a collection of these 
3 oscillators. Therefore, it was surprising to see that when the parameters are tuned as 
mentioned above, two of the masses coalesce and one gets a higher derivative oscillator in 
addition to a normal one. The model at that limit is at most Jordan-diagonalizable. It 
is a ghost-free higher-derivative model whose properties were put under scrutiny recently 
[9, 10, 12], where quantization was also carried out. 

In this work, we have only considered the non-interacting theory. It would be quite 
interesting to consider the lowest order interacting theory to see whether the higher derivative 
model stays ghost-free (or as far as the classical solutions are concerned, whether there are 
stable regions or not). Such an analysis was carried out by Smilga [11] in a simplified model. 



APPENDIX A: THE PAIS-UHLENBECK OSCILLATOR 



Here we would like to recapitulate the nonrelativistic limit of the higher derivative field 
theories considered a long time ago by [8]. Pais and Uhlenbeck studied a generic scalar field 
theory model whose Lagrangian is of the form 



£ = - 



1 



N 



1=1 



Here we will take N = 2 and examine only the nonrelativistic limit, which basically means 
the dropping down of the spatial derivatives. Thus consider a Lagrangian of the form 



r -'- f ••2 / 2 I 2\ -2 I 2 2 2I 
2^ ~*^'^l+'^2)9 +^1^2? I 



(Al) 



and note that uji 7^ UJ2 and Ui = 002 cases differ drastically. We will call the latter case as 
the degenerate Pais-Uhlenbeck oscillator which is similar to the case studied in this paper. 
The naive Ostrogradski Hamiltonian simply reads 



1 



- 222 



(A2) 



which gives the following, not necessarily positive, energies for solutions of the equation of 
motion 

UJi 7^ U2 



E 



al'^i - ^2 



^M + 4)-cjI{cI + cI) 

L0{C5C6 - CyCg) +4 + 0? 



(A3) 



CUi = UJ2 



where Cj(i = 1, . . . , 8) are constants. To be able to quantize these theories, a proper Hamil- 
tonian in terms of momenta can be defined by using Dirac's constraint analysis. Taking 
q — 'K-i one gets 



H 



1 



2 2\ 2 

,^1 + '^2)^ 



(A4) 
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where and Pg are the usual canonical momenta. One can see the glimpses of two coupled 
oscillators here, but one still needs to "diagonalize" this Hamiltonian. For Ui = uj2, just like 
the oscillator studied in the bulk of this paper, the Hamiltonian is only Jordan diagonalizable. 
For oji UJ2: one can use the following Pais-Uhlenbeck variables 

Qi = q+\, Q2 = q+\. (A5) 
In this case the Hamiltonian becomes 

H = \^^M + ^IQI) - \^^M + ^IQI) , (A6) 

which is a collection of two oscillators except that one has the wrong sign. For a discussion 
of the quantization of the lox ^ LO2 case, we refer the reader to [9, 10]. 



APPENDIX B: THE SOLUTIONS 

Let us now examine the steady-state solutions of the D — Z Pais-Uhlenbeck oscillator. 
For this purpose, we assume that the mass parameter m has been scaled appropriately and 
the field (or Z2 thereof) has been redefined accordingly such that the problem reduces to 
investigating the solutions of the (V^— m^)^<^ = equation on the (2-spatial dimensional full) 
plane. Below we will give the exact solution of this equation in terms of highly complicated 
special functions. However, one can easily extract the asymptotic behavior of the solution 
at spatial infinity, which is simply (/? ~ (ci + C2r)e~'"''. Of course, we have assumed "nice" 
boundary conditions and discarded the divergent parts to find this. 

First consider the simpler equation (V^ — m'^)(p = 0. For the problem at hand, one 
can assume circular symmetry, take (p — (f{r,9), employ a 'separation of variables' ansatz, 
introduce a new independent variable x = mr, arrive at the well known modified Bessel 
equation and, further assuming the solutions to be free of any angular dependency, write 
down 

^{r,9) = R{mr) = R{x) = aIo{x) + (3 Kq{x) , 

for some integration constants a and (3. Thus, our original problem reduces to examining 
the solutions of the following ordinary differential equation 

(PR 1 dR „/ X ^ ^, / X 

if we further demand that the solutions of (V^ — m?Yip = are finite and regular at large 
X (or r) . One can view this as a two-point boundary condition problem, where the solution 
R{x) is required to be finite and bounded for both small and large values of x. Using 
standard Green's functions' techniques, one arrives at 

(PX PCX) \ 

Koix)Ioiy)Koiy)ydy + j h{x) Kl{y) y dy\ , 

which can be evaluated using the software package Mathematica. The outcome of this is 

R{x) = ip{mr) = -(5 {^-x" h{x){-Kl{x) + Kl{x)) + ^Ko{x)u{x) 
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where 



u{x) = MeijerG[{{3/2}, { }}, {{1, 1}, {0}}}, ^1 - ^ / ' r(^^-l)^^^^' ^"'^ ' 



in which the contour of integration on the s-plane is between the poles of r(— s — 1/2) and 
r(s + 1). By using the properties of the F-functions, one can also write 



and note that T{—s — 1/2) has poles at s = —1/2, 1/2, 3/2, . . . and r(s + 1) has poles at 
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